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Abstract 

Using a novel numerical spectral method, we have constructed an AdSg-CFT4 solution to the Einstein equation with 
a negative cosmological constant A that is asymptotically conformal to the Schwarzschild metric. This method is 
independent of the Ricci-DeTurck-flow method used by Figueras, Lucietti, and Wiseman. We have perturbed the 
solution to get large static black hole solutions to the Randall-Sundrum II (RSII) braneworld model. Our solution 
agrees closely with that of Figueras et al. and also allows us to deduce the new results that to first order in 1 /(-AM 2 ), 
the Hawking temperature and entropy of an RSII static black hole have the same values as the Schwarzschild metric 
with the same mass, but the horizon area is increased by about 4.7/(-A). 
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1. Introdcution 

The Randall-Sundrum II (RSII) braneworld model 
is one of the braneworld models suggested for solv- 
ing the hierarchy problem. The braneworld model in- 
cludes a higher-dimensional spacetime that is called 
the bulk and a lower-dimensional spacetime called the 
brane, which is embedded in the bulk. All matter and 
fields in the standard model are supposed to propagate 
on the subspace manifold, the brane, but gravity is the 
only force that can propagate through the whole space, 
the bulk. The RSII model is a warped five-dimensional 
braneworld model with an extra dimension that can be 
large. It is a very important question whether station- 
ary large black holes exist within the RSII braneworld 
model. If not, then the extremely strong observational 
evidence for very-nearly stationary astrophysical black 
holes would be a nearly conclusive reason for reject- 
ing that theory as possibly physically realistic. Af- 
ter various conjectures and claims J2I |3[ B El 01 that 
large black holes do not exist in the Randall-Sundrum II 
(RSII) braneworld model JJ], Figueras and Wiseman (01 
(henceforth FW) recently found such solutions by per- 
turbing an AdS5-CFT4 solution that Figueras, Lucietti, 
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and Wiseman ]8[] (henceforth FLW) had found earlier by 
Ricci-DeTurck flow. This AdS5-CFT4 metric is a solu- 
tion to the Einstein equation with a negative cosmolog- 
ical constant A that is asymptotically conformal to the 
Schwarzschild metric. Because the Schwarzschild met- 
ric appears at an AdS^ boundary with an infinite scale 
factor, it may be viewed as a black hole of infinite mass. 

We had independently searched for and found the 
infinite-mass black hole solution by a different numeri- 
cal method and were preparing to perturb it to get large- 
mass RSII black hole solutions when the Figueras et al. 
papers appeared. Here we report that our numerical so- 
lution agrees well with that of Figueras et al. and thus 
adds further evidence for the existence of large RSII 
black holes, despite the doubts expressed by previous 
work. 

We used a spectral method, expressing the compo- 
nents of the 5 -dimensional metric in terms of Legen- 
dre polynomials in the two nontrivial coordinates, with 
the appropriate boundary conditions imposed. Then we 
chose the 210 coefficients of the polynomials to mini- 
mize the integrated square of the error of the Einstein 
equation, finding that we could reduce this by eight 
orders of magnitude from the case with no free pa- 
rameters (constant polynomials). The integrated square 
is based on the Gauss-Legendre quadrature numerical 
method, and the minimization procedure uses the sim- 
plex search method for multivariable functions. This 
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strongly suggests that we are numerically near an ex- 
act solution, though of course our limited computational 
resources meant that we could not use an infinite num- 
ber of parameters to reduce the numerical error all the 
way to zero. This approach to solving Einstein equa- 
tions is novel, and the good agreement of our results 
with the Figueras et al. results illustrates the success of 
the method, especially in comparison with the failure of 
various previous numerical attempts. 

We present an explicit approximate metric for the 
black hole on the brane. Using this approximate met- 
ric we demonstrate that the area of an RSII black hole 
on the brane is slightly greater than a black hole in pure 
four-dimensional general relativity, and to leading or- 
der, the relations between the mass, Hawking tempera- 
ture, and Bekenstein-Hawking entropy are precisely the 
same as in four-dimensional general relativity. In other 
words, astrophysical-sized black holes in the Randall- 
Sundrum II braneworld model are extremely close to 
what four-dimensional general relativity would predict. 
Although this may be disappointing for those hoping 
to distinguish between the two models by observations 
of black holes, it it highly encouraging for those who 
have postulated that the Randall-Sundrum II braneworld 
model is a model consistent with our observations. 

2. Infinite Black Hole Metric 

For brevity of notation, we use units in which the 5- 
dimensional bulk cosmological constant is A — -6. We 
start with the AdS-black string metric 

ds 2 = -^-r \dw 2 + U(ry l dr - U(r)dt 2 

+ r 2 d£l 2 \ (1) 

where U(r) = 1 -2M/r and where dSl 2 = dfi+sin 2 6d(f> 2 
is the unit two-sphere metric. Letting r = 2M/y and 
w = 2M/v gives 
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The hypersurfaces of constant v are Schwarzschild 
metrics of mass m(v) = v/2. The curvature at y > 
diverges at v = 0, so this black string metric is singular. 
We modify the metric by adding some y 2 terms to re- 
move this singularity, and we also introduce four metric 
functions to give 



ds 1 = A 



dv 2 



v 2 + y 2 



+ B 



(v 2 +y 2 )dy 2 
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- 4C(v 2 + y 2 ){\ - y)dt 2 + D—dQ 2 . (3) 
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We then replace v, which ranges from to <x>, by x = 
y 2 /(y 2 + v 2 ), so that the metric becomes 



ds 2 = A(1-jc) 
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where < x < 1, < y < 1 and A, B, C and D are 
smooth functions of x and y. The coordinate bound- 
aries are these: x = is the infinite AdS boundary that 
is conformal to the Schwarzschild metric when we im- 
pose A-B-C-D-l there, y — is the ex- 
tremal Poincare horizon, x — 1 is the axis of symmetry 
where the two-sphere shrinks to zero size and where we 
impose A - D for regularity, and y — 1 is the black 
hole horizon where we impose the regularity require- 
ment B — C. 

The most general metric satisfying all the symme- 
tries for our problem has five components. On the other 
hand, since the metric functions depend non-trivially on 
the two coordinates, x and y, and the choice of these 
is gauge dependent, one can reduce the number of the 
metric components to three. The common method for 
finding a unique solution for the Einstein equation nu- 
merically is to fix the gauge before discretization. Oth- 
erwise, one will have a family of solutions parametrized 
by one function that is gauge dependent. But in our 
case, we assume four unknown functions instead of 
three as mentioned, and we still get a unique solution. 
Our explanation for this result can be related to our re- 
striction of A(x, y), B(x, y), C(x, y), and D(x, y) to poly- 
nomials for simplicity; we also tried rational functions, 
but they did not seem to work numerically so well. Hav- 
ing polynomials of some fixed finite order means that 
with such restricted functions, some gauges are better 
than others. For clarification, we can consider the case 
of a spherically symmetric static metric 



ds 2 = -A(x)dt 2 + -^—dx 2 + -J—dtl 2 . 
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If we consider the restriction for A(x,y), B(x,y), and 
C(x,y) to be polynomials, then with having all three 
functions, one can find A(x) = 1 - x, B{x) — x 4 - x 5 , 
and C(x) = x 2 solves the vacuum Einstein equations. 
However, if one choose the gauge B(x) = 1, no polyno- 
mials of finite order would give an exact solution, and 
we would expect greater error. On the other hand, we 
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are not looking for an exact solution, so with our restric- 
tion to have a fixed order of polynomials for each func- 
tion, surely we would get a better result with more func- 
tions, even if for an exact solution one or more functions 
would be just gauge. 

We impose these regularity conditions and also solve 
the Einstein equation to lowest order in x by writing 

A = \-x{\-x){\+2f{y))+x 2 g{y) 

+ x 2 (l -x)A{x,y), 

B = \+xf(y)+x 2 B{x,y), 

C = l+xf(y)+x 2 B(x,y)+x 2 (l-y)C(x,y), 



D = \+x{\-x)(\+f(y))+x 2 g(y) 
+ x 2 (l-x)D(x,y). 



(6) 



With units such that A = -6, the vacuum Einstein 
equation in the 5-dimensional bulk is 



E a f3 = R a p + 4g aj3 - 0. 



(7) 



We define the integrated square error of the Einstein 
equation to be 



< 5 W 5 x, 



(8) 



where we choose At = 2n in order to get a definite fi- 
nite integral (assuming that E a pE al3 falls off fast enough 
toward the infinite AdS boundary at x = 0, where the 
metric determinant (5, g oc l/x 6 diverges, so that the in- 
tegral converges). 

We choose polynomials for the functions f(y), g(y), 
A(x,y), B(x,y), C(x,y), and D(x,y) and numerically 
vary the coefficients to minimize the integrated square 
error /. For A = B = C= D = 1, I * 4038, but 
when we went up to sixth-order polynomials with a total 
of 210 coefficients, the integrated squared error was re- 
duced to 0.00004238, nearly eight orders of magnitude 
smaller. The maximum value of the squared error at 
any point within the 5-dimensional spacetime was then 
E a pE a P = 0.000154. Thus we appear to have strong ev- 
idence that our numerical method is converging toward 
an exact solution of the infinite black hole metric. 

Because our metric uses different coordinates from 
those used by FLW, it is not easy to make many com- 
parisons over the entire bulk 5-dimensional manifold. 
We have found that the minimum value for the length 
scale given by the inverse fourth root of the total trace of 
the square of the Weyl tensor, {C a p y sC al3yS y 1 ^ ', is about 
0.206 in our metric, which is within 4% of the value 
0. 198 that FLW privately reported to us from their met- 
ric. However, we shall make many more comparisons 
below for the 4-dimensional large black hole metric. 



3. Finite Black Hole Metrics 

To get a black hole metric with a large but finite mass, 
we need to replace the infinite AdS boundary at x = 0, 
where the metric is conformally Schwarzschild but with 
infinite mass, by an RSII brane with induced metric y^ v , 
and with a mirror image of the bulk metric on the oppo- 
site side of the brane. Assuming no matter on the brane, 
in our units with A = -6, the Israel junction condition 
is [Kfjy] = -2y MV , where the square brackets denote the 
difference in the second fundamental form from the 
back side to the front side, where it has the opposite 
value. Hence K^ v = -y^ on the front side, where we 
are using the opposite sign convention from FW. 

We write the asymptotic form of the bulk metric near 
the infinite AdS boundary as 



ds 2 = -^[dz 2 + g llv (z,x)dx' 1 dx v j, 



(9) 



where the z is the exponential of the negative of the out- 
ward proper distance as one approaches the AdS bound- 
ary at infinite proper distance (where z — > 0), and where 
x denotes the other four coordinates (not just the single 
x coordinate used above). Then we use just enough of 
the Fefferman-Graham (FG) expansion jfTol |TTl llSl Ivh 



to write 



~g, v (z,x) * gfXx) + z 2 [-U^(x) 



+ z L v (x). 



(10) 



Here gjj(x) is the conformal metric on the infinite AdS 
boundary, and f jUV (x) is a tensor that is covariantly con- 
served in the conformal metric gffi(x) and which gives 
the second set of constants of integration for the bulk 
Einstein equation when written as second-order equa- 
tions in z. We extracted polynomial forms for t MV (x) 
at the boundary from our sixth-order polynomial met- 
ric coefficients for the bulk metric. 

The infinite-mass black hole bulk solution is equiv- 
alent to putting the RSII brane at z = and setting 
gfvi x ) - 8% h ' tne Schwarzschild metric with unit hori- 
zon radius. Since the Schwarzschild metric is Ricci-flat, 
Rfy(x) vanishes, so then there is no z 2 term in the FG ex- 
pansion. One can also show that then t MV (x) is traceless 

as well as conserved en HUES]. An approximation 

for it can be extracted from the numerical data for the 
bulk solution found above. 

For a large but finite mass black hole, we put the 
brane at small z — e and perturb the conformal met- 
ric at z = to g^y = + e 2 /J Alv . This gives the Ricci 
tensor a small perturbation, so that now it is of order e 2 , 
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and the z 2 term in the FG expansion is no longer zero. 
The perturbation in g^ also perturbs t^x) away from 
its traceless form, but since that term is multiplied by z 4 
in the FG expansion, for results to lowest nontrivial or- 
der in e, it is sufficient to use the original value of t MV (x) 
from the infinite-mass black hole bulk solution. 

Using the fact that = (z/2)d z [g,j V (z,x)/z 2 ], to 
lowest nontrivial order in e the Israel junction condi- 
tion implies that R®>(x) - (l/6)R^g^(x) * 4e\ v (x). 
Since to this order t MV (x) is traceless, one further gets 
that the Ricci tensor of the perturbed conformal metric 
8 { °v(x) = gl? + e 2 V is B$(x) * 4e\ v (x), which with 
a knowledge of t MV (x) is sufficient to determine h^y and 
hence the spherically symmetric static metric gfy(x). 
Then the induced metric on the brane is 



1 



1 



(11) 



The bulk Einstein equation plus the Israel junction 
condition for a brane with the RSII value of the ten- 
sion and without matter imply that the Ricci scalar of 
the brane metric is zero. We can achieve this to first 
order in the perturbation for a generic static spherically 
symmetric metric on the brane by going to the gauge 
h' t — (which is equivalent to choosing the coordinate 
y so that on the brane g t , = -(1 - y) after rescaling 
t so that g tt = -1 at radial infinity, y = 0), defining 
h e „ = fit = (y 2 /6)F(y), and then setting 



,v 2y\l-y) { T7 , dF 

«« = ~z~t~. — ^— \F+y— 



3(4-3?) 



dy 



(12) 



If we now define (2M) 2 = 1/e 2 = 6/(-Ae 2 ) after 
reverting to general units in this last expression, and if 
we define a new radial coordinate p = 2M/y, then to first 
order in e 2 = (3/2)/(-AM 2 ), we can write the metric 
on the brane (after rescaling the time coordinate t by a 
factor of AM) as 



l ds 2 = y^ydx^dx 
1 



1- 



p-2M I dF 
{-Ap 2 )p-l.5M\ ~ P ~dp 
2M\ 
P I 



l-fjdp 2 



dt 2 + 



P 2 + 



1 



(-A) 



dQ. 1 



(13) 



One can show that the asymptotic behaviour of t MV (x) 
(which goes as 1/p 5 for p » 2M with known co- 
efficients JH 0]) implies that F approaches unity as 
p — > oo or y — > 0. To fit the FLW numerical data 
f^J; w (jc) which they kindly sent us, and to fit 



(2) 

our numerical data ,(x) 
and F2 = Fpuw to be cubic polynomials in y = 2M/p 



f™'(x), we took Fi = F om - 



with the constant coefficient set to unity and then chose 
the other three coefficients in each case to minimize the 
respective 



Ji = 



Jp'A^yp^d'x 



(14) 



where for each of the two values of i (i = 1 for the FLW 



data and i — 2 for our data) in the numerator A/S = 



f™ - tfy is the difference between the t^l(x) given by 
the cubic for F,(x) and the t^l given by the numerical 
data. The integral in the denominator was included to 
make J a normalized mean-square error and has tj^ w {x) 
given by the FLW numerical data, the same in each case 
to give a constant normalizing factor. The factor of p 4 
was included to increase the weight of the large-p part, 
though the integrals are still dominated by the small-p 
part, since t^ v {x) drops off asymptotically as the inverse 
fifth power of the radial coordinate p la LZD . 

For the FLW numerical data f^[r w (x) (which was con- 
strained to be traceless and proved to be very nearly 
conserved and matching the predicted y 5 dependence at 
small y), J\ was minimized at 7flw ~ 0.0000620 for 



FLW 



1 _ 1 .062^U .554^f 



\ P I 



P I 



0.120 

P 



(15) 



For our numerical data f™'(x) (which was not quite 
traceless and conserved and also had a small spurious 
y 4 term), the normalized mean-square error J2 was min- 
imized at 7 our 0.00 139 for 



(2 f^) +0 .434(^ 



0.059 



2M 
P 



(16) 



From the fact that J om « 22 7flw, clearly our data is 
less accurate than the FLW data, which is not surpris- 
ing since we varied only 210 parameters in our spectral 
method, whereas FLW used grids of 40 x 40 (or 1 600 
points) and of 160 x 160 (or 25 600 points). Also, the in- 
dividual coefficients of these two cubics have large rel- 
ative differences, but the ratio of the two cubics them- 
selves never differs by more than 1.3% from unity, so 
they show good agreement between what is generated 
by our numerical data and by what is given by the FLW 
data. 

We also used the integral of Eq. (TBI with At^y = 
t^y"' - f^Jr w (x), the difference between the stress tensor 
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2M 2.5M 3M 4M 6M 10M25Moo 
P 

Figure 1: At the top is an llth-order polynomial fit F\\ to the FLW 
data, which gave normalized mean-square error J\\ = 0.0000572, 
92% of /flw- Because the diff erences from F\\ of the cubic fits Fflw 
and F ow are too small to show up when plotted directly on this graph, 
at the bottom we have expanded these differences by a factor of 50 
and plotted 50(Fflw _ F\ 1 ) (bottom curve) and 50(F O ur — ^ll) (middle 
curve). 



& m generated by our F our fit to our data and the stress 
tensor tj^ w (x) given directly by the FLW data. This 
gave J 3 = /our fit vs. FLW data ~ 0.000214 m 3.4 7 FLW , so 
the t^° m (x) generated by our F our (x) fits the FLW data 
about 6.5 times better than it fits our data. This is also 
not surprising, since the fy 1 " (x) generated by our F 0UI (x) 
fit was constrained to be both traceless and conserved, 
whereas the f°" r (x) extracted directly from our data was 
not. 

Furthermore, we calculated the mean-square error be- 
tween the t^y"'(x) generated by our F 0U i fit to our data 
and the f^ LW (jc) generated by the F FLW fit to the FLW 
data, using in Eq. <Q3) Af^ = t^T - t^ w (x), and got 
J4 = Jam- fit vs. flw fit * 0.000146 ~ 2.4 7 FLW , so the 
tav m (x) generated by our F om - fits the t^ LV (x) generated 
by the Fflw fit to the FLW data nearly 9 times better 
than it fits the t°" r (x) directly extracted from our data, 
which is not quite traceless and conserved, as the f^™' (x) 
generated by the fitting F om - is constrained to be. 

We also calculated the ratios between the values of 
each of the three individual components of tj^ w (x) as 
given by the FLW data, the fit to the FLW data given by 
^tv LW W> an d the fit to our data given by fy m (x). These 
ratios were generally within 1-2% of unity, with the 
maximum differing by less than 2.9%. The ratio of the 



hy's generated by F our and by Fflw, which involves a 
derivative of F as given in Eq. (11), did differ by up 
to about 9.3%, so when one is neither near the black 
hole horizon (p = 2M exactly in our gauge) nor near 
2M/p = (at both of which limits = 0), the devi- 
ation of g pp from the Schwarzschild value may not be 
given to very high precision by our F om . However, the 
other results appear to agree within a very few percent 
from those given by FLW, which gives strong indepen- 
dent confirmation of their results. 

One can readily calculate from the metric (12) for any 
F that has only a constant term and negative powers of 
p, as Fflw given by Eq. (14) and F om given by Eq. (15) 
do, that the ADM mass is precisely M and that the sur- 
face gravity of the black hole horizon is exactly 1 /(AM), 
the same as for the Schwarzschild metric. Of course, 
aside from the numerical approximations for determin- 
ing F, the metric (12) is only correct to first order in 
our perturbation parameter 1 /(-AM 2 ), so there might 
be corrections to the surface gravity of a static RSII 
black hole to second order in 1 /(-AM 2 ). However, one 
can deduce that to first order in 1 /(-AM 2 ), the Hawking 
temperature and entropy for the RSII black hole have 
the same values as they do for the Schwarzschild met- 
ric. 

On the other hand, the horizon area is shifted from 
the Schwarzschild value As c h = An(2M) 2 to Arsh = 
An[(2M) 2 + F(l)/(-A)], where F(l) is the value of F 
on the horizon, at y = 2M/p = 1. The fit to the 
FLW numerical data gives fpLw(l) ~ 0.372, and the 
fit to our numerical data gives F our (l) » 0.373, which 
agree within about 0.3%. Therefore, the change from 
a Schwarzschild black hole to an RSII black hole on a 
brane with the same ADM mass M increases the hori- 
zon area (but not the Hawking entropy) by the amount 
AA = 47rF(l)/(-A) « 4.67/(-A), where here we used 
the FLW data value as probably more accurate. (The 
value from our data would give a coefficient of about 
4.69 rather than 4.67.) 

4. Conclusion 

We have provided independent numerical evidence 
in support of the numerical discovery of Figueras and 
Wiseman |7[] of the existence of large static black holes 
in the Randall-Sundrum II braneworld model JH, by 
a significantly different numerical method. Our re- 
sults agree quite well with theirs, such as giving an in- 
crease in the black hole horizon area over that of the 
Schwarzschild metric of the same mass by 4.69/(-A) 
or 4.67/(-A) respectively. (This is a new result, not re- 
ported in |0].) We have obtained a good closed-form 
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approximation to the metric of the black hole on the 
brane, Eqs. (12) and either (14) or (15). We have also 
shown the new result that to first order in our pertur- 
bation parameter 1 /(-AM 2 ), the Hawking temperature 
and entropy of the black hole is the same as that of a 
Schwarzschild black hole of the same ADM mass M. 

If large black holes did not exist in the Randall- 
Sundrum II braneworld model, the astrophysical obser- 
vations of such black holes would have been strong ev- 
idence against the viability of that model. However, 
our confirmation of the large black holes in RSII found 
by Figueras and Wiseman J3l, and the fact that they 
are very nearly the same as Schwarzschild black holes, 
show that the RSII model is not excluded in this respect. 
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Figueras, James Lucietti, and Toby Wiseman and 
greatly appreciate their sharing their detailed numeri- 
cal data with us for comparison. CC acknowledges an 
Avadh Bhatia Postdoctoral Fellowship at the University 
of Alberta. This research was also supported in part by 
the Natural Sciences and Engineering Research Council 
of Canada. 
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